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Abstract 

The scalar, vector, and tensor components of the (generalized) deuteron 
electric polarizability are calculated, as well as their logarithmic modifications. 
Several of these quantities arise in the treatment of the nuclear corrections to 
the deuterium Lamb shift and the deuterium hyperfine structure. A variety of 
second-generation potential models are used and a (subjective) error is assigned 
to the calculations. The zero-range approximation is used to analyze a subset of 
the results, and a simple relativistic version of this approximation is developed. 
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1 Introduction 


The spectra of hydrogenic atoms and ions have been measured to such high precision 
that nuclear properties play a significant role [1]. This is especially true of the deu¬ 
terium atom, which has large nuclear contributions to both the Lamb shift[2, 3] and 
the S-state hyperfine structure[4], These large nuclear contributions result from the 
weak binding of the deuteron and from the concomitant large size of this simplest of 
all nuclei. 

The Lamb shift can be defined as the difference between calculated energy level 
values and those from predefined reference values, which is typically a Dirac spectrum 
modified by reduced-mass effects. This shift results after removing hyperfine split¬ 
tings, which is equivalent to having a spinless nucleus. Hyperfine structure results 
from the electron’s interaction with the nuclear spin, and can be linear in that spin 
(the usual type) or quadratic (quadrupole hyperfine structure). This typically means 
that the effective nuclear interaction with the electron can be a scalar, a vector, or a 
tensor in character. 

The dominant electron-nucleus interactions (beyond the point-nucleus Coulomb 
potential) occur in first-Born approximation in their electromagnetic coupling. Exam¬ 
ples are the scalar (or L=0) nuclear finite-size modification of the Coulomb potential, 
which dominates nuclear effects in the Lamb shift, and the tensor (or L=2) modifica¬ 
tion that dominates the quadrupole hyperfine structure. The leading-order (vector) 
hyperfine structure results from the electron’s interaction with the nuclear magnetic 
moment. 

Sub-leading nuclear effects are generated by two-photon processes, which necessar¬ 
ily involve a loop integral over a virtual photon momentum. The Lamb shift integral 
has a weaker dependence on that momentum (i.e., the process is “softer”) than the 
corresponding hyperfine process (which is therefore “harder”). Because the most 
important parts of the two-photon nuclear amplitudes involve sequential electron- 
nucleus electromagnetic interactions, numerically important contributions arise from 
inelastic intermediate nuclear states, and thus are non-static or polarization phenom¬ 
ena. These mechanisms can also be scalar, vector, or tensor in type. 

The leading-order nuclear contribution to the Lamb shift is proportional to the 
mean-square nuclear charge radius, while the subleading order is proportional to the 
nuclear electric polarizability and its logarithmic modification. The polarizability is 
defined in terms of the nuclear electric-dipole operator, D, and the fine-structure 
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constant, a, as 


2a |(jV|Q|0>| 2 

E 3 ko E " ~ E« ’ 

and its logarithmic mean-excitation energy [2], E, by 

log(2 E/m e ) = log [2(E n - E 0 )/m e 

o - -no 


( 1 ) 

( 2 ) 


Both are purely non-static, and involve the (virtual) excitation of negative-parity 
intermediate states with energy, Ejy, from the deuteron ground state with energy, E 0 . 
The precise form of the argument of the logarithm has been dictated by conventional 
atomic physics usage. 


The dominant part of the sub-leading-order nuclear (vector) hyperfine structure 
is determined by “Low” moments[4], while smaller subleading-order contributions 
involve ordinary dipole excitations of a different type. Adopting a uniform convention 
for constant factors that will be apparent later we define 

= (0|D|Af) x (Af|D|0) = -Ao\D x 5|0>, (3) 

6 N^O 6 

and 

2a- 

log(2£/m e ) a = -i— £ (0\D\N) x (N\D |0) log [2 (E N - E 0 )/m e ] . (4) 

6 N^O 

The factor of i is necessary to make the appropriate part of a real. We also note 
that in non-rclativistic approximation D x D vanishes because all components of 
D commute with each other. This is not true when relativity is incorporated in 
the calculation (see Eqn. (34) of [5], which demonstrates that D x D is imaginary). 
In that event the commutator contributes an essential part of the derivation of the 
Drcll-Hearn-Gerasimov sum rule[6] and the Low-Energy Theorem[7]. The quantity a 
is more appropriately described as a sum rule, rather than as a polarizability. 


The leading-order contribution to the quadrupole hyperfine structure is driven 
by the quadrupole part of the ordinary Coulomb interaction, whose scale is set by 
the nuclear quadrupole moment, Q. A sub-leading contribution is determined by the 
tensor (L=2) component of the nuclear electric polarizability tensor 


^ _ 2a (0\De\N)(N\D a \0) 

O'E n / x r-i r-i 

o Nf L 0 — -n o 


(5a) 


which has been calculated only once long ago using separable potentials and the RSC 
potential model[8]. Those calculations were used in a purely nuclear (as opposed to 
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molecular) determination[9] of the deuteron’s quadrupole moment. No estimate of 
the polarizability effect on the determination of the deuteron’s quadrupole moment 
from the HD molecular hyperhne structure has ever been made. There will also be a 
logarithmic modification of this tensor polarizability, which is obtained by inserting 
log [2 (E N — E 0 )/m e ] in the sum over intermediate states. 

A system with spin 1 can in general have scalar, vector, and tensor polarizabilities, 
and these can be of the ordinary type defined in Eqn. (5a), or of a type without the 
factor of (En — Eq) in the denominator, determined by 

D fio = ^-(0\D^D a \0 ), (5b) 

which was illustrated by Eqn. (3) (and more accurately called a sum rule). Each 
of these 6 types can have a logarithmic modification. Because the deuteron is an 
important nucleus and deuterium an important atom, we will calculate all twelve 
of these polarizabilities (or sum rules) using modern second-generation potentials. 
Only three of the polarizabilities have been calculated before, and only two with 
modern potentials. None of the logarithmic modifications of the sum rule above 
have been previously calculated. Our approach will be non-relativistic in keeping 
with the potential models we use. Our calculations will be based upon Podolsky’s 
method[10] for computing second-order perturbation-theory matrix elements, and on 
several integration tricks, one of which has been used in the past. 


2 Tensor Polarizabilities 


The electric polarizability is most easily calculated using second-order perturbation 
theory and the coupling of the nuclear electric-dipole-moment operator to a uniform 
electric field, which leads immediately to Eqn. (5a). This equation is fully equivalent 
to 

< 4 ° = y <0|B' 3 |A*“}, (6a) 


where [10] 


(H — E 0 ) |AT“) = D a |0) 


( 66 ) 


is solved subject to hnite boundary conditions. Note that D does not connect the 
ground state (the only bound state) of the deuteron to itself. Resolution of Eqns. (6) 
into partial waves is necessary in order to perform a numerical calculation. 


Because we wish to compute all tensor components of the electric polarizability 
tensor we present a very brief derivation of the partial-wave decomposition that up- 
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dates older work[11, 2], We begin with Eqns. (5a) and (6) and convert them to a 
scalar by contracting the Cartesian index a with a constant vector E a and the index 
/3 with a constant vector E' 13 . The deuteron initial state dC = | SM) (and analogously 
for the final state) depends on the azimuthal quantum number M and we formally 
remove it from the problem by dehning[12] a (vector) projection operator (appropri¬ 
ate for S = 1): e\ x {M) = S\ t M- Performing the sum l-S'A) results formally 

in a scalar quantity independent of spin projections and amenable to manipulation. 
A similar projector J2y e iy(M') (S\'\ is used for the final (i.e., leftmost) deuteron 
state. 


We define general orthonormal spin-angular wave functions for the deuteron sys¬ 
tem 

<l>jM=( Y e(r)®Xi)jM, (7) 

which couples the usual angular wave function that depends on the direction r of the 
internucleon vector f to the (unit-) spin wave function, y. The deuteron’s full wave 
function is then given by 

I’d = ^ a t( r ) <l>\ i(0> ( 8 ) 

e=o,2 

where cto(r) = u{r)/r and a, 2 (r) = w(r)/r expresses ai in terms of the conventional 
radial wave function components of the deuteron. We have suppressed for now the 
azimuthal quantum number, M. 


We ignore the tiny effect of the n-p mass difference on the deuteron center-of-mass 
(CM), and find D = |(rf — rf )r*expressed in terms of the isospin operators r* of 
nucleon i. The dipole isospin operator generates T = 1 states when acting on the 
deuteron, together with a residual numerical factor of | that we will ignore until 
later. If the right-hand-side (RHS) of Eqn. (6b) (contracted with E a and e* x (M)) is 
expanded in terms of (f/j functions, one Ends 


where 


RHS 


I] <!> L j-{E l ®e 1 )j 

L= 1,3 
J=0,l,2 




9j( r ) 

V3 


(9) 


Qq =u — \[2w (10a) 

g\ =u + w/V2 (10 b) 

g-l — u — v^2 w /10 (10c) 

gl =3v / 3w/5 (lOd) 


are the relevant radial functions. We manipulate the left-hand-side of Eqn. (6b) 


5 



(similarly contracted) into the same form as Eqn. (9) 

i AS ,« )B « = z^ £ in) 

r L= 1,3 V3 

J= 0,1,2 

where fi is the n-p reduced mass and the functions fj(r) satisfy 

(H L J - E d )^-f}(r) = gi . (12) 

We note that for total angular momentum J = 2 the L — 1 and L = 3 orbital 
components are coupled by the tensor force. We do treat that coupling properly, 
although it is not reflected in the simplified notation employed in Eqn. (12). 


The matrix element in Eqn. (6a) (including its factor of (2a/3), a factor of 1 from 
each of the two dipole operators, the factors of (—2 fi) and l/\/3 from Eqn. (11), and 
the l/y/3 from Eqn. (9)) then becomes 



X (£> e-;v (£i ®«,)., (-i)'E 4, 

J L 


(13a) 


where 


CL j — 


-fid 

~9~ 


drr 2 fj{r) g L j{r) 


(136) 


now expresses the entire content of the electric polarizability tensor in terms of pro¬ 
jection operators and matrix elements. Equation (13a) is not a convenient form, and 
we recouple it so that the projectors of the same type are coupled together. We also 
note that only the J = 2 part of the a/j terms has two non-vanishing components. 
We therefore define 


and 

which yields 


where the (real) quantities 


— a j + 5j t 2 a \, 

(E[ ® Ei) k ■ {£.'* ® ei) K , 

a E a -»■ 


(14) 

(15) 

(16) 


6° — g(^0 + 3Al + 5 A 2 ) 

b\ = —(— 2Aq — 3Ai + 5 A 2 ) 
6 

62 = — 3^4, + A 2 ) 

6 
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(17a) 

(176) 

(17c) 



determine the tensor properties of the nuclear physics. 


The structure of the X K operators corresponds to tensors of order k in both the 
Cartesian indexes a and /3, and in the effective (azimuthal) spin dependence. The 
dependence on the two spin-projection operators e\ and ej* is indeed equivalent to 
using the Wigner-Eckart Theorem [12] on the nuclear matrix elements. This allows 
us to rewrite the spin factors \ K as effective operators in the nucleus total-angular- 
momentum Hilbert space, determined by powers of the angular-momentum operator, 
S of the nucleus: 


13a 

Oi E = OiE 


£ + fae ^ + T 


'S a S l3 + S f} S a 25 a/3 ' 


(18) 


The three coefficients a E , cr, and r are, respectively, the scalar, vector, and tensor 
components of the polarizability. Equation (1) (a trace) defines a E , while cr, and r 
are defined by 



eV3 {SS\D*\N){N\D»\SS) 

N^O ~ ^0 


(19a) 


and 


r = 3a| 3 - a E , 


(196) 


where the deuteron should be in the state M’ — M — S — 1. Finally the relationships 


a E = b 0 (20a) 

a = — b\ (206) 

r = - b 2 (20c) 

determine the various polarizabilities in terms of the partial waves. Note that cr 
is real, since the Cartesian vector cross product differs from the spherical one by a 
factor of —i\f2 and the Wigner-Eckart Theorem guarantees that the spherical result 
is overall real. 


The sum-rule quantity D 3a is decomposed in strict analogy to Eqn. (18), with 
coefficients, s, v, and t replacing a E , a, and r. In the non-relativistic approximation 
that we employ (or, equivalently in the deuteron, the impulse approximation) the 
former quantities are related to conventional deuteron moments by ^ = f (r 2 ) c , 
v = 0, and t = ^-Q, where (r 2 ) c / l is the mean-square charge radius and Q is the 
quadrupole moment. Note that in the representation of Eqn. (18) the quantity a in 
Eqn. (3) becomes a = v S. 
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3 Logarithmic Sum Rules 


Calculations of the logarithmic modification of the basic polarizabilities or sum rules 
use the trick of adding an arbitrary energy £/ to the energy denominator in Eqn. (5a), 
where / has the dimensions of energy and £ is dimensionless. This defines 


_ 2a (0|£>^|Ar)(jV|£>°|0) 
E 3 fe, (f + E N - E 0 


( 21 ) 


We first integrate a§“(£) with respect to £ from 0 to A, where A is very large compared 
to any relevant energies ( E N — E 0 ). This produces 

/ /V<4“(0 = ^ y(0|B' J |lV)(iV|B“|0)log(A//(B A , -£„)]. (22) 

6 N^O 


We split the integration region from [0, A] into [0,1] plus [1, A], In the second region 
we change variables from £ to l/£. We also note that .///a d£/£ = log A. Putting 
everything together we find that A can be taken to infinity if we use 


^ E {0\D l3 \N)(N\D a \0} log [2(E n - E 0 )/m e \ = log (2E/m e ) = 

d 0 


- I #/<*?(«)- /‘f 

Jo £ 


~D ri “ 


+ D j3a log (2//ro e ), (23) 


where D 1301 is dehned in Eqn. (5b). 


A similar set of manipulations was developed previously in which the integral 



2a (0|L> /3 |AT)(iV|T) a |0) 


o / -/ 

d TV 


E]\f — En 


log [(E n — Eq)/ ef) 


(24) 


was split into the regions [e, 1] plus [l,oo], and the integration variable for the second 
region was also changed to l/£. This led to a special case of 

of(0) log (2 E/m.)= [' ^[a£ , (fl-ag’(0)+af , (l/fl]-af*(0) log(m«/2/). (25) 

Equations (23) and (25) provide a tractable scheme for calculating logarithmic mod¬ 
ifications of our basic polarizability, a^“(0), and of the sum-rule quantity D ,3a . The 
results are independent of the scale parameter, /. Note also that since the various 
polarizabilities are particular linear combinations of partial waves, the logarithmic 
modifications are also the same linear combinations involving those partial waves, 



and this is a convenient way to perform the calculations. We will also see in the 
next section that the zero-range approximation provides an excellent starting point 
for understanding the scalar polarizabilities (the vector and tensor ones are signifi¬ 
cantly smaller and vanish in this approximation). In the Appendix we develop the 
zero-range form of Eqn. (5a), which allows the analytic calculation of the logarith¬ 
mic modifications to a e and D aa . We also develop a zero-range model based on the 
Relativistic Schrodinger Equation (RSE). 


4 Results and Conclusions 


We have calculated the scalar, vector, and tensor components of both and D 3a , 
together with their logarithmic modifications. These calculations were performed with 
seven different second-generation potential models, including the Argonne Vis (AV18) 
[13], the Reid Soft Core 1993 (RSC93), and 5 Nijmegen models[14, 15], including the 
full model (no partial-wave expansion) as well as the local and non-local Reid-like 
models. The last two types had been fitted to both relativistic and non-relativistic 
forms of the deuteron binding energy. 

It has been known for a long time that the deuteron mean-square radius and 
the electric polarizability are rather accurately predicted by the zero-range approxi¬ 
mation (see the Appendix). That approximation over-predicts the polarizability by 
approximately 1% and D aa by less than 2%, and is therefore an excellent starting 
point for investigating the uncertainties in the four scalar quantities. The largest 
uncertainty in the zero-range results is due to As, the asymptotic S-wave normal¬ 
ization constant, whose value was determined in phase-shift analyses[16] to be As = 
0.8845(8) frn -1 / 2 . This leads via Eqn. (A6) to the zero-range result, af = 0.6378(12) 
fm 3 , and via Eqn. (A2) to _D" a = 0.01916(4) frn 2 , where we use the relativistic form 
of the deuteron binding energy for both quantities. 

The second-generation potentials are sufficiently accurate that they can be re¬ 
garded as alternative phase-shift analyses. We therefore expect that the values of the 
electric polarizability and D aa will scatter around a central value with the variance 
of 2 parts/thousand associated with A|. In order to check whether the potential- 
model variation of the remaining 1% of a# and roughly 2% of D aa is small, we have 
scaled each calculated quantity by (A]g xp /A™ ode1 ) 2 , and have examined the remaining 
variations. The scaled values are listed in Table I. This procedure verihes that the 
dominant uncertainty in the scalar quantities is the error in As- 

The tensor quantities should be expected to scale like rj, the asymptotic D/S 
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Table I 


Table 1: Scalar, vector, and tensor components of D ,3a (viz., s, v, t ) and (viz., «e, 
a, r), followed by the product with the appropriate logarithmic factor, log(2i7/m e ), 
for each case. Results were calculated using a number of different potential models as 
discussed in the text, and the “error bar” results from combining a subjective estimate 
of the spread in the results after scaling to the experimental values of A? s and rj, as 
discussed in the text, with the variances of those two quantities. All calculations used 
the impulse approximation for the dipole operator and assumed equal-mass nucleons. 
Note that the first two rows have been multiplied by a factor of 10 to make the entries 
more uniform in size. 


type 

Scalar 

Vector 

Tensor 


D fio x 10 

0.1882(4) 

0.00000 

0.01322(10) 

fm 2 

D^ a log {2E/m e ) x 10 

0.6327(12) 

0.0003(1) 

0.0503(4) 

fm 2 

/3a 

a E 

0.6330(13) 

-0.00092(5) 

0.0317(3) 

fm 3 

c4 Q l°g (2 E/m e ) 

1.8750(36) 

-0.0023(2) 

0.1014(8) 

fm 3 


(amplitude) ratio of the deuteron. Those quantities in Table I have therefore also 
been scaled by ^ ex P/^ model ; where the experimental value of rj was determined from 
phase-shift analyses[16] to be rj exp = 0.0253(2). This scaling considerably reduces the 
scatter in the calculated tensor results. The small size of i] also roughly accounts for 
the size of the tensor quantities relative to the scalar ones. The result in Table I is 
nearly 10% smaller than the one in Ref. [8], reflecting a smaller modern value of // 
[16], 

The vector quantities are very suppressed and are sensitive to details of the nuclear 
force that vary from model to model. They are sensitive to the forces in the 3 Pj states, 
and in the absence of those forces can be shown to be determined by the square of the 
D-state wave function. The quoted uncertainties are inferred solely from the model 
variations. Note that the vector part of D !3a vanishes identically in non-relativistic 
approximation, and this accounts for the very small sizes of the vector quantities. 

In summary, we have calculated a variety of polarizabilities and sum rules for 
the deuteron that are generated by (unretarded) electric dipole interactions. These 
quantities have been divided into scalar, vector, and tensor components, and include 
logarithmic modifications of each. Our numerical techniques allow us to generate 
all such components at no extra cost. The scalar polarizabilities were previously 
calculated[2], and have not changed significantly. The vector components of the 
sum rule play a role in the ordinary deuterium hyperhne structure[4] and are very 
small. The tensor components determine part of the deuterium quadrupole hyperhne 
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structure [1], although calculations of this effect have not yet been performed. The 
sizes and uncertainties of the various quantities were analyzed using the zero-range 
approximation and various scales appropriate to the deuteron. 


5 Appendix 


The zero-range approximation is motivated by the asymptotic dominance of radial 
matrix elements that contain (positive) powers of the distance between the nucleons. 
The mean-square radius is an obvious example, as is the electric polarizability, since 
each is weighted by two powers of the inter-nucleon separation. The usual version of 
this approximation is to assume that in intermediate states the nucleons lie outside 
the range of the nuclear force (i.e., we set that force to zero) and in the initial and 
final deuteron states we ignore the D-wave and use the asymptotic form of the S-wave 
function. The zero-range deuteron wave functions are therefore given by 


u ZI (r) = A s exp (-nr ), (Ala) 

w zr (r) — 0 , (Alb) 


where Ag is the asymptotic S-wave normalization constant, k = y/2/iEb, /i is the 
n-p reduced mass, and E b is the (positive) deuteron binding energy. Ignoring the 
tiny difference in the proton and neutron masses (the neutron lies slightly closer to 
the CM than the proton) one finds that the mean-square charge radius and ( D 2 ) are 
proportional in impulse approximation (which we have assumed as a consequence of 
the non-relativistic approximation), but not otherwise. Recalling the factor of (2a/3) 
from Eqn. (5b) and two factors of | from the two dipole operators, the zero-range 
value of D aa is given by 


D 


aa 

zr 


® A S 

24 k 3 


0.01916(4) fm 2 , 


(A2) 


where the numerical result uses the experimental value of Ag and the relativistic value 
of k defined below. The calculation of the remaining quantities requires an analytic 
expression for the scalar electric polarizability with the £/ insertion. It is convenient 
to choose / to be the deuteron binding energy (i.e., Eb = —Eq), so that in momentum 
space we have 


£E b + E n — Eq 


(1 + Q k 2 + f 

2 n 


(A3) 


and the required Green’s function is therefore a simple modification of the usual 
zero-range Green’s function for the deuteron: 


Go(0 


2/i 

47rr 


exp (— nr ) 


(A4a) 
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k = \Jl + £ k . (A4b) 

Performing the integrals in Eqn. (5a) using Eqn. (Al) for the wave function and 
Eqn. (A4) for the Green’s function leads to the zero-range result 

afK ) = aMi (* 2 + + 4 kk) 

E ^’ 12k 3 (k + k) 4 V 7 

This gives the well-known result 

of (°) = - 0.6378(12) fm 3 , (A6) 

where the numerical result uses the experimental value of As and the relativistic 
value of k dehned below. Equations (22) and (24) can now be used to calculate the 
logarithmic modifications of a^ and ( D 2 ) zr . The logarithmic modification of a^ is 
determined by 

log (E/.E b ) = log 4 — ^ , (A7a) 

or -A = 2.23214- • while the modification of (D 2 ) zr is determined by 

log {E/Eb) = log 4 — ^ , (A7b) 

o 

or -p- = 3.38592 ■ ■ -. Both results are very simple and quite accurate. 


Our final task is to estimate the size of one class of relativistic corrections to a 2 ^. 
We begin with the so-called Relativistic Schrodinger Equation, which we construct 
for two non-interacting nucleons with identical masses, M, by summing the kinetic 
energies of each: 


E 


(M 2 + Pi) 2 - (M 2 + p 2 ) 


H 

2 \ 2 


T = 0 


(A8) 


A potential could also be added to the kinetic energy. In the center-of-mass frame of 
the two particles (with momenta p and — p, respectively) 


2 (M 2 + p 2 ) 


T = 0, 


(A9a) 


indicating that the energy of a bound deuteron would be given by 

l 

E d = 2 (M 2 - k 2 ) 2 = 2 M - E b . (A9b) 

Since the rightmost (experimental) result holds in all cases, it clearly makes a dif¬ 
ference if the non-relativistic approximation 2 M — K 2 r /M is substituted for the 
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square root in Eqn. (A9b). For this reason we have labelled the relativistic value 
of k as K r and the non-relativistic approximation as K nr . They are related by 
K r = K nr (l ~ rf lr /8M 2 ). 

1 

Equation (A9a) also holds if we multiply it by E cm +E' cm , where E' cm = 2 [m 2 + p 2 ^j 2 . 
This reduces that equation to non-relativistic form , but with K r replacing K nr . Equa¬ 
tions (Al) therefore still hold mutatis mutandis. This does not apply to the Green’s 
function, however. If we invert Eqn. (A9a) and multiply top and bottom by the iden¬ 
tical factor E cm + E' cm , the denominator has the desired form —4 (k 2 + p 2 ), but the nu¬ 
merator now contains the factor (E cm + E' cm ), which we rewrite as 2E cm + (E' cm — E cm ) 
and note that the second part of this expression cancels a similar factor in the denomi¬ 
nator. The remainder is very short ranged (range ~ 1 /M) when Fourier transformed, 
and in keeping with the zero-range approximation we ignore this term. Thus the 
appropriate zero-range Green’s function for the RSE is simply (1 — k^/M 2 ) 1 ^ 2 Go, 
where the form in Eqn. (A4a) holds if we replace 2/i by M and k by K r . Thus pre¬ 
vious results for o?e hold is we use K r everywhere for k and multiply by the factor of 
(1 — k 2 /M 2 )^ 2 . This produces a correction (compared to 1) ~ — k 2 /2M 2 ~ —0.0012, 
which reflects the expected size of relativistic effects for the denteron. For a full 
treatment, see Ref. [17]. 
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